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Abstract 


Throughout the world; topographic maps are compiled by manually operated 
stereoplotters that recte&te the geometry of two wide-angle overlapping 
stereo frame photographs. Continuous imaging systems such as strip 
cameras, electro-optical scanners, or linear arrays of detectors (push 
brooms) can also create stereo coverage from which, in theory, topography 
can be compiled. However, the instability of an aircraft in the atmosphere 
makes this approach impractical. The benign environment of space permits 
a satellite to orbit the Earth with very high stability as long as 
no local perturbing forces are involved. Solid-state linear-array ; -.nsors 
have no moving parts and create no perturbing force on the satellite. 
Digital data from highly stabilized stereo linear arrays are amenable 
to simplified processing to produce both plauimetric imagery and elevation 
data. A satellite, called Mapsat, including this concept has been 
proposed to accomplish automated mapping in near real time. Image 
maps as large as 1:50,000 scale with contours as close as 20-m interval 
nay be produced from Mapsat data. 

Background 


The geometry of stereo mapping photographs, whether taken from aircraft 
or satell’.te, is well known and documented. Transforming such photographs 
into topographic maps is a relatively slow and expensive process that 
for many critical steps defies automation. Compared to an aircraft, a 
satellite offers the unique au vantages of much greater stability and 
uniform velocity. 

Utilizing these advantages, a sensing system in space can now provide 
imagery of mapping quality, even Chough a continuous electro-optical 
imaging system is used instead of a mapping camera with its inherent 
high geometric fidelity. The next gen.^ration of space sensors will 
include solid-state linear arrays (fig. 1) that involve no moving 
parts. By continuous Imaging with very high geometric fidelity they 
will permit, at least in part, the automated mapping of the Earth 
from space in three as well as two dimensions. The fundamental difference 
between conventional ana continuous stereo methods is illustrated 
by figure 2. 


* Aoproved for publication by Director, U.S.G.S. 
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At least four papers have been published that relate directly to autouated 
three-dimensional mapping. In 1952, Kaez (1) showed how height 
measurements could be made with a stereoscopic continuous-strip camera. 

The geometry of such a strip camera and stereo Hr ear arrays Is basically 
the same. In 1962, Elms (2) elaborated on the strip camera concept 
and Indicated Its advantages over frame cameras as a possible component 
of an automated mapping system. In 1972, Helava and Chapelle (3) described 
the developtaenc of Instrumentation by which a conventional stereomodel 
can be scanned using the eplpolar-plane* principle, and thus reducing 
image correlation from a two-dimensional to a basically one-dimensional 
task. 

In 1976 Scarauo and" ffrumm' ’described the autdnatexl stereo-mapper 
AS-llB-X which utilizes the eplpolar-scan concept and one-dimensional 
digital Image correlation described by Helava and Chapelle. Thus the 
concept of reducing photograametric data stero correlation iron o 
to one dimension Is well established. The cited literature, however 
does not describe the possibility of imaging the Earth directly in 
stereoscopic digital form suitable for one -dimensional processing. 

Beginning in 1977 a serious effort to define a stereo satellite cr 
Suereosat (5) was vndertaken by NASA. The Stereosat concept calls for 
linear-array sensors, looking fore, vertical and aft, but Its principal 
objective Is to provide a stereoscopic view of the Earth rather than 
to map It in automated mode. There ...a other ways of obtaining stereo 
Imagery with linear arrays. The French SPOT (6) satellite can look 
left or right of the track and thus achieves stereo by combining Imagery 
from nearby passes of the the satellite. NASA's Multispectral Linear 
Array (MLA) concept (7), as so far defined, calls foi fore and aft 
looks through the same set of optics by use of a rotating mirror. 

However, neither the SPOT nor NASA's MLA approach are considered 
optimum for stereo mapping of the Faith, as neither Is designed to 
acquire data In continuous fort.. 

Mapsat Geometric Concept 

Linear arrays represent a relatively new remote sensing concept. Five 
papers on this subject were presented at the ASP/ACSM a.inual convention 
during March 1978 (8,9,10,11,12). These papers concentrated on detector 


*An eplpolar plane is defined oy two air or space exposure (imaging) 
stations and one point on the ground. 
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technology and the application of linear array sensors in a vertical 
imaging mode. W^'sh (13) rece:.tj.y aescrlbed the geometry of linear arrays 
in stereo modv, although his error analysis for such a system is based on 
measuxements made from images rather than computations based on the digital 
data* 

3y combining che technology/ of linear arrays, the concept of epipolar- 
plane scanning, and the experiv ace gained from Landsat and other space 
sensing systems, Mapsat was defined (14), and its proposed parameters are 
listed In Table 1. The Mapsat concept was the work of several individuals, 
but perhaps the single most important contribution was that of Donald Light 
(verbal communic*;''ion) , then of the Defense Mapping Agency, who first 
suggested that epipolar planes, as described by Helaya (3) and used in the 
AS-113^X plotter, could he achieved directly from space and that topographic 
data might then be extracted in real time. There are several feasible 
configurations by which linear array sensors can continuously acquire stereo 
data. It vas decided tnac the system must permit selection from the three 
spectral oands, provide for two base-to-height ratios of 0.5 and 1.0 and 
be compatible with the epipola. concept. Figure 3 illustrates the 
configuration selected to accomplish the stereoscopic as well as monoscopic 
functions. 

Acquiring stereo data of the Earth in epipolar form directly from space is 
the fundamental geometric concept of Mapsat. The epip >lar conditions 
shown in Figure 4 implies that five poxnts — the observed ground point ?, 
the two exposure stations S and S . and the two image detectors f and a — lie 
in a single plane. 1 this epipolar condition is maintained as the satellite 
moves along its orbit, every point P observed by detector f in the forward 
looking array will also be observed subsequently by detector a in the aft 
looking array. Thus image correlation can be obtained by matching 
the data stream from detector f with that from a — a one-dimensional 
correlation scheme. This description applies equally to the use 
of the vertical with either the fore-or aft-looking array but involves 
a weaker (0.5) base-to-height ra.io than the described use of the fore 
and al'" arrays (base-to-height ratio of 1.0). In practice the data streams 
from more chan one detector tcay be involved *ince there will normally be some 
offset in the oath of a given pair of detectors. Moreover under certain 
conditions, cotielacion may be improve' by a limited expansion of 
the correlation function to two aiDens'ov.3. 
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Becau«e each detector array is looking at a different portion of the 
Earth at any given time, Earth rotation complicates the epipoLar 
condition. As shown in figure 5, this complication can be overcome 
by controlling the spacecraft attitude. This description is obviously 
simplified; further complications involve such factors as the ellipsoidal 
shape of the Earth, variations in the orbit, spacecraft stability, 
and even very large elevation differences. The spacecraft position and 
attitude- must be precisely determined by such- systems as the Global 
Positioning Sv^cem (GPS or NAVSTAR) and frequent stellar referencing. 
Satellite attl.ude control involves gyros and inertial wheels, and, 
when a satellite is free of perturbing forces created by moving (actuated) 
parts, attitude can be aiaintalnad for reasonable periods to the arc-second. 

Of course, the sensing system nust retain precise geometric relationship 
to Che attitude control system. Defining the correct satellite attitude 
and the rates in yaw, pitch, and roll to maintain the epipolar condition 
requires precise raatheraarical analysis. T%’o independent analyses, one 
by Kowell of ITEiC (15) and the other by Snyder (16) of U.S. 

Geological Survey, confirm Mapsat's geometric feasibility, and a 
U.S. patent has been allowed on the concept. Table 2 Indicates the 
maximum deviations from the epipolar condition caused by the varions 
expected error sources. This table is based on a half orbit (5U minutes) 
which covers the daylight portion to which imagery is basically Limited. 
Attitude rate errors would be considerable if only corrected once evep-* 

50 minutes but, as Che table indicates. iC-minute intervals based on 
stailar reference reduce the errors to a reasonable amount. Ten-minute 
stellar referencing using scar sensors as described by Junkins et al. , 

(17), is conslderei reasonable. Computer programs have been developed 
that result in the epipolar plane condition being maintained as long as 
adeo”.ate positional and attitude reference data are available and proper Iv 
utilized. Figure 6 illustrates the simplicity of elevation decerraination 
in an epipolar plane which is the key element of Mapsat. 

Obviously, the Mapsac concept can be effectively implemented only if 
stringent specifications regarding orbit, stability, reference, and sensor 
systems a ' met. Table 3 lists the MapsaC geometric requirements as define 
to date, and each is considered to be within the state of the arc. 

Mapping Accuracy 

By meeting the geometric requirements indicated and achieving stereo 
correlation, the resulting tna,. accuracy is compatible with scales as large 
as 1:50,000 and contours as -lose as 20 m interval based on U.S. National 
Map Accuracy Standards. Reference 15 covers this analysis in some detail. 
Such accuracies result from the indicated geometric requirements and the 
following factors: 
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o Linear e^ray dececCocs are poslcioned wlch sub'nticroa accuracy* 

o Optical distortion effects', when accounted for by calibration, 
are negligible. 

o .Atmospheric refraction, because of the steep look angles, is of 
a very low order and i« reasonably well known; air-to-water 
refraction is also known where underwater depth detemination 
is Involved. 

0 Relative timing, which is referenced to data acquisition, is 
accurate to within the microsecond. 

0 Digital stereo correlation-,' where uniquely achieved, provides 
three dlnenslohal ’rbot-mean-*square (rms) positional accuracy 
to within half the pixel dimension. 

These considerations result in relative positional errors for defined 
points of only 6 to 7m (rtos) both horizontally and vertically. This 
vertical accuracy requires the 1.0 base-to*-helght ratio. Such accuracy is 
adequate for the mapping indicated but assumes chat control is available 
for reference to the Earth's figure. As indicated by ITEK (13) and the 
author (19), control points of 1,000 km spacing along on orbital path will 
be adequate for such a purpose. Where no control exists the absolute accuracy 
of the resultant maps, with respect to the Earth's figure, may be In rms 
error by 50 to 100 m although their internal (relative) accuracy remains 
at Che 6 Co 7 m rms level. 


Stereocorrelation 


The determination of elevations from stereo data requires the 
correlation of Che spectral -esponse from the sane poini; or group of points 
as recorded from two different positions. In the aerial photography case 
these two positions are the camera stations, whereas with linear arravs in 
space the two recording positions are constantly moving with the satellite. 

In the photography case, correlation is achieved bv orierting the two 
photographs to model the acquisition geometry. Once this is done, correlation 
can be achieved by the human operator, or the image stereoraodel can be scanned 
ant correlated by automateu comparison of the signal patterns from the two 
photographs. A system such as the AS-llB-X (3,4) generates cne-dimens tonal 
digital data in epipolar planes from the model, ''n theory, eptpolar data ^ 
should be correlated much faster than that from a system that must 
search In two dimensions to establish correlation. In practice, the automated 
correlation oi digital dat^ has been only partially successful; and, as Mahoney 
(18) has recently pointed out, correlation by either manual or automated systems 
is still a slow and costly process. To date, no one has acquired original sensor 
data in epipolar form. -s, no one can really say how well such data can be 
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auCooatlcally correlaCed, until a satellite such as Hapsat is flown. 

Siaulation using digitized aerial photographs or linear-array stereo-sensing 
of a terrain model are relevant experiments worth conducting. However, they 
will provide only partial answers, since the degree of correlation will depend 
on the area involved. The characteristics of the Earth's surface, coupled 
with related conditions, such as the atmosphere and Sun angles, are highly 
varied; which means that the degree of correlation will also be highly varied. 
This problem does not imply that the Mapsat concept has not been validated. 
Having stereo data organized in linear digital form is of obvious advantage 
to create the three-dicsnsional model of the Earth's surface. !lany areas 
will correlate in one-dicensional node, others will require two-dimensional 
treatment, and still other areas nay not correlate at all. By properly, 
defining the satellite parameters and data processing,, the correlation 
f unctloc can he" opetalzed aad ralsed welL above that obcainable- from 
wide-angle photography systems. For example, digital data can readily be 
modulated to enhance contrast or edges that make up the patterns on which 
correlation depends. Photography can also be nodulated, but it is far more 
difficult vand less effective) chan digital-data modulation, as fi.la lacks 
the dynamic ranga and sensitivity of solid-state detectors. Mapsat will 
acquire data in an optimum form for automated correlation, which will expedite 
the precise determination of elevations and create digital elevation data 
that are becoming a basic tool for many disciplines. 

Acquisition Modes and Products 

As previously described <'.4), Mapsat is designed to be operated in a 
wide variety of nodes. These include variation in resolution (lO-n elements 
on up), spectral bands, swath width, and stereo modes. Such flexibility 
permits optimum data acquisition without exceeding a specified data- 
transmission rate that is now defined at 48 megabits per second (Mb.-'s). 

The Earth's surface is highly varied, and data product requirements 
are likewise highly varied. By varying the acquisition modes and, in turn, 
producing a variety of products, the data management problem beconec 
complicated as compared to existing systems such as Landsat which produces 
only two bai.lc types of data. However, solving this data management problem 
is a small price to pay for a system that can meet a wide variety of 
requirements for remotely sensed data of the Earth. Only four primary 
products are expected from Mapsat as follows: 

Ca) Raw-data digital tapes from which quick-lbok images can be 
displayed in near real time. 

(b) Processed digital image tapes calibrated both radiometri tally 

and geometrically to a defined map projection. Such data will be 
two-dimensional ( planimetrlc) but describe the Earth's radiance 
(brightness) in multispectral form as is now accomplished by 
Landsat Multispectral Scanner tapes. 
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(c) Processed digital tapes, again calibrated both radiometrlcally 
and geometrically, but which now describe the Earth’s surface 

in three dimensions (topographically) with an asscciated radiance 
value* Such tapes are, in effect, digital elevation data sets of 
the Earth's surface. 

(d) Standardized images, both black*and-white and in color, which 
include geometric corrections and radiometric enhancements* Such 
corrections and enhancements will be of recognized general value 
and of a type that can be performed without undue delay or 
excessive cost* The images would also be of standardized scale* 

Scow these: four haeic produeta, a: .wide, variety of derivatives can be made whiah 
include the following: 

(a) Black-and-white and multicolor image maps and mosaics at scales 

as large as 1:50,000, or even 1:25,000 (1:24,0(0) where nap accuracy 
standards are not required* 

(b) Thematic displays and maps involving such subjects as land cover 
and land use classification. 

(c) Maps which depict the Earth's topography by such means as contours 
(as close as 20-m interval), slopes, elevation zones, shaded relief, 
and perspective display. 

Ccnciusion 


Mapsat will not meet all anticipated remote sensing requirements, 
and it will in no way replace those air-photo surveys required to meet 
mapping requirements for scales Larger chan 1:50,000 and contour intervals 
of less than 20 ra* What it will do, is provide a precise three-dimensional 
mulcispectral model of the Earth at reasonable resolution and in digital 
form. Moreover, the satellite will record the changing responses of the 
Earth's surface as long as it is In operation. 

Mapsat can be built today at what is considered be a reasonable 
cost (15) as It is based on available components and technology/. Moreover, 
it is d?^signed for simplified operation and data processing- Assuming that 
a*' operational Earth-sensing system will be flown, surely Mapsat is a 
deserving candidate for such a job* 
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Mapsat Parameters 


0 Orblc — Same as Landsat 1, 2 and 3 (919 km alc^. 

o Sensor — Linear Arrays~—Three optics looking 23° forward, vertical 
and 23° aft* Three spectral bands: 

blue green 0.17 — 0.57 um 
red 0.57 - 0.70 um 

near 18. 0.76 ^ 1.05 um 

o Swath — 130 km or portion thereof. 

o Resolution““Variable~Dawn to 10 m element. 

o Transmi;isioo — S (or X) band, compatible with Landsat receivers 
but with rates up to 13 Mb/s. 

o Processing-One dimensional. Including stereo. 


Table 1 
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Mapaat Geometric Requireaf-t. *s 

o Positional Determlnatloa of Satelllte~lO to 20 in 
all cnree axes. 

2 / 

o , Pointlag Accuracy—tflthln— 0.1 of vertical. 

2/ 

o Pointing Determination — Within— 5 to 10 arc seconds 

o Stabtlicy of Satellite— ^tatiohai rates'Vithiai'^ 10’^ 
degrees /second. 


IJ rms (Icr ) 

very high probability (3a) 


Table 3 
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DETECTORS 



Figure 





Lintar arrays of datactora In 

fowl IG to 30 cn long 




HapsAC SansoT Coaflsarmcl^ (aot ta sc2l«^ . 

Opclcj A, a. and 0 4T« * rigid part of satcUita. 
Optic 3 sanaca tha sama atrip 60 sacoada aftar A; 
optic C, UO sacoada aftar A. Aay coobinacioa of 
A, 3, and C producaa scarao. Optica A and C art 
of abc\it ICZ loagtr focal laagtii to provida 
raaolution compaciPit vit^ optic 3. 
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figure 
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